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The short and intermediate range order of an amorphous Ge3oSe7o alloy produced by Mechan- 
ical Alloying were studied by Reverse Monte Carlo simulations of its x-ray total structure factor, 
Raman scattering and differential scanning calorimetry. The simulations were used to compute the 
Gg^g e (r), G%f_g{r) and G^g{r) partial distribution functions and the Sg c M § c (^), S^g(K) and 
Sft^SciK) partial structure factors. We calculated the coordination numbers and interatomic dis- 
tances for the first and second neighbors and the bond-angle distribution functions Oiji(cos8). The 
data obtained indicate that the structure of the alloy has important differences when compared to 
alloys prepared by other techniques. There are a high number of Se-Se pairs in the first shell, and 
some of the tetrahedral units formed seemed to be connected by Se-Se bridges. 

PACS numbers: 61.10.Eq, 61.43.Bn, 05.10.Ln, 87.64.Je 
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I. INTRODUCTION 

Amorphous semiconductor materials like chalcogenide 
glasses present a great potential for application in tech- 
nological devices, such as optical fibers, memory materi- 
als and switching devices, but their use is limited due to 
several factors. One of them is the difficulty in obtain- 
ing information about atomic structures, which define 
the short-range order (SRO) of the alloy. In this context, 
the structures of amorphous Ge^Sei-^ (a-Ge^Sei-^) and 
liquid Ge^Sei-a; (?-Ge x Sei_ x ), in particular Ge33Se67 
(GeSe2), have been extensively studied by several exper- 
imental techniques, like neutron diffraction (ND)i*Sii4, 
x-ray diffraction (XRD)£, extended x-ray absorption fine 
structure (EXAFS)£*£ and Raman spectroscopy (RS)&&. 
On the theoretical side, molecular dynamics simulations 
(MP) 10 ! 11 ! 12 ! 13 : 14 ! 15 : 16 have been carried out to under- 
stand the SRO in these liquids and glasses in terms of 
two possible and distinct models. In the first one the 
distribution of bonds in the structure is purely randomic 
and determined by the local coordination numbers and 
composition. In the second one, there is a strong SRO 
and the structure is formed by well defined structural 
units, e.g., corner-sharing GeSe 4 / 2 (CS) tetrahedral and 
edge-sharing Ge 2 Se 8 / 2 (ES) bitetrahedral units. The dis- 
tribution of these units gives raise to a medium, or inter- 
mediate, range order (IRO), whose signature is the ap- 
pearance of a first sharp diffraction peak (FSDP) in the 
neutroniSia or x-ray structure factors^ at many composi- 
tions. In particular, ND experiments performed on melt- 
quenchedM (MQ) GeSe 2 (MQ-GeSe 2 ) showed a FSDP in 
the total structure factor S(K) which was associated with 
correlations in the range of 5-6 A. As described in Ref.Q, 
this alloy is formed by CS and ES units with heteropolar 
bonds but there are homopolar bonds in very small quan- 
tities. These results are reproduced by MD simulations in 
Z-GeSe^ 12 : 13 ! 15 and a-GeSe^ 10 i n i 14 i 16 , except the FSDP in 



the Bathia-Thorton (BT)il concentration-concentration 
structure factor Scc(K). It should be noted, however, 
that almost all available data about a-Ge^Sei-a; alloys 
were determined for MQ samples, and the preparation 
method can affect the SRO and IRO, as it was reported 
by Takeuchi et al& by comparing the structures of a- 
Ge3oSe7o produced by MQ and by vacuum evaporation 
(VE) techniques. Tani et aim- studied the a-GeSe 2 pro- 
duced by Mechanical Grinding (MG) of its crystalline 
counterpart and also found some structural differences. 
These differences are important because some physico- 
chemical properties can be altered and improved as de- 
sired by choosing an appropriate method of preparation. 

It is well known that the Mechanical Alloying tech- 
nique (MA)i£ introduces a high quantity of disorder and 
defects in the structure of the materials produced us- 
ing this process. Thus, in this paper, we investigated 
the SRO and IRO of an amorphous Ge3oSe7o alloy pro- 
duced by MA (MA-a-Ge3oSe7o) starting from the elemen- 
tal powders of Ge and Se using Raman spectroscopy, Dif- 
ferential Scanning Calorimetry (DSC), X-ray diffraction 
and reverse Monte Carlo (RMC) simulationa 19 ! 20 : 21 ' 22 of 
its XRD S(K ). We were interested in studying two main 
points. Since it is not obvious that starting from the 
powders of Ge and Se and submitting them to a milling 
process we would obtain an alloy formed by (ordered) 
structural units like CS or ES units, first of all we would 
like to know if the alloy produced by MA contains these 
units. Besides that, even if these units are formed the 
SRO and the IRO of the alloy can be significantly altered 
by the high quantity of defects and disorder introduced 
by the MA process when compared to MQ samples, for 
instance. Therefore, the second point is to determine the 
local structure of the alloy itself, finding coordination 
numbers and interatomic distances. At our knowledge, 
this is the first time that such study is reported concern- 
ing an a-Ge3oSe7o alloy produced by MA. 
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II. THEORETICAL BACKGROUND 
A. Structure Factors 



sion terms and Ci is the concentration of atoms of type i 
(A is the radiation wavelength). The partial reduced dis- 
tribution functions dj (r) are related to Sij (K) through 



1. Faber and Ziman structure factors 

According to Faber and ZimanS 3 ., the total structure 
factor S{K) is obtained from the scattered intensity per 
atom I a (K) through 



S(K) = 



I a (K) [(f(K)) {f{K)f 



n 71 
i=l j=l 

where K is the transferred momentum, Sij(K) are the 
partial structure factors and Wij (K) are given by 



2 f°° 

G lj (r) = - / K [S^ (K) - 1] sm(Kr) dK . 
11 Jo 



From the Gy (r) functions the partial radial distribution 
function RDF^- (r) can be calculated by 



RDFy(r) = 4 7 rp c 3 r 2 + rG tJ (r) . 

where po is the density of the alloy (in atoms/A 3 ). Inter- 
atomic distances are obtained from the maxima of Gij (r) 
and coordination numbers are calculated by integrating 
the peaks of RDFy(r). 



and 



Wij{K) 



cjCjUK)fj(K) 



i 



Here, f l {K) = f 0i (K) + //(A) (A) is the atomic scat- 
tering factor, //(A) and /"(A) are the anomalous disper- 



2. Bathia and Thornton structure factors 

The BT structure factors can be related to the FZ 
onesii. For a binary alloy the BT number-number struc- 
ture factor <Snn {K) is given by 



S NN (K) = clS n (K) + 4S 22 (K) + 2 Cl c 2 S 12 (K) , (1) 
where Sij(K) are the FZ partial structure factors and 
Ci is the concentration of element i. The BT number- 
concentration structure factor iSnc(^0 i s 



Svc(K) = c lC2 { Cl [Su(K) - S 12 {K)] - c 2 [S 22 (K) - S 12 (K)} } 



(2) 



r 



and the BT concentration-concentration structure factor Sec [K) is found through 
I 



Sec (K) = ci c 2 { 1 + a c 2 [S n {K) + S 22 (K) - 2S 12 (K)] } 



(3) 



r 



RMC Method 



The basic idea and the algorithm of the standard 
RMC method are described elsewherai2*22i2i*22i and 
its application to different materials is reported in 
the literaturo 24 ! 25 ! 26 ! 27 . 28 ! 29 ! 30 ! 31 ! 32 ! 33 ! 34 ! 35 ! 36 ! 37 ! 38 . In the 



RMC procedure, a three-dimensional arrangement of 
atoms with the same density and chemical composition 
of the alloy is placed into a cell (usually cubic) with pe- 
riodic boundary conditions and the G^ M (r) functions 
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corresponding to it are directly calculated through 

n RMC (r) 
J A~k par z Ar 

where n^ MC (r) is the number of atoms at a distance 
between r and r + Ar from the central atom, averaged 
over all atoms. By allowing the atoms to move (one at 
each time) inside the cell, the G^ MC (r) functions can be 
changed and, as a consequence, S RMC (K) and S nMC (K) 

are changed. Thus, S KMC (K) is compared to the S(K) 
factor in order to minimize the differences between them. 
The function to be minimized is 

v> 2 = i£[s(^-s RMC (^)] 2 , ( 4 ) 

i=l 

where the sum is over m experimental points and S is re- 
lated to the experimental error in S(K ). If the movement 
decreases ip 2 , it is always accepted. If it increases ip 2 , it 
is accepted with a probability given by exp(-AV> 2 / 2 ); 
otherwise it is rejected. As this process is iterated ip 2 de- 
creases until it reaches an equilibrium value. Thus, the 
atomic configuration corresponding to equilibrium should 
be consistent with the experimental total structure fac- 
tor within the experimental error. By using the G^ M (r) 
functions the coordination numbers and interatomic dis- 
tances can be calculated. In addition, the bond-angle 
distributions &iji(cos9) can also be determined. 

III. EXPERIMENTAL PROCEDURE 

The MA-a-Ge3oSe7o alloy was produced by consider- 
ing a binary mixture of high-purity elemental powders 
of germanium (Alfa Aesar 99.999% purity, particle size 
< 150 /jm) and selenium (Alfa Aesar 99.999% purity, 
particle size < 150 /im) that was sealed together with 
several steel balls into a cylindrical steel vial under an 
argon atmosphere. The ball-to-powder weight ratio was 
5:1. A Spex Mixer/Mill model 8000 was used to perform 
MA at room temperature. The mixture was continuously 
milled for 40 h. A ventilation system was used to keep 
the vial temperature close to room temperature. The 
XRD pattern was recorded in a powder Siemens diffrac- 
tometer equipped with a graphite monochromator, using 
the CuK Q line (A = 1.5418 A). The total structure factor 
S(K) was computed from the XRD pattern after cor- 
rections for polarization, absorption, and inelastic scat- 
tering, following the procedure described by Wagner—. 
The /' and /" values were taken from a table compiled 
by Sasaki 4 ^. Raman measurements were performed with 
a T64000 Jobin-Yvon triple monochromator coupled to 
a cooled CCD detector and a conventional photon count- 
ing system. The 5145 A line of an argon ion laser was 
used as exciting light, always in backscattering geometry. 



The output power of the laser was kept at about 200 mW 
to avoid overheating the samples. All Raman measure- 
ments were performed at room temperature. DSC mea- 
surements were taken with a heating rate of 10°C/min 
using a TA 2010 DSC cell in a flowing argon atmosphere. 



IV. RESULTS AND DISCUSSION 
A. Raman Scattering 

FigureQlshows the RS spectra of MA-a-Ge3oSe7o, c-Se 
and c-Ge. The alloy has bands at around 195, 215, 237 
and 255 cm" 1 . The band at 195 cm -1 is associated with 
the Ai breathing mode of CS units and the 215 cm" 1 
band is related to the A\ breathinglike motions (com- 
panion peak) of Se in the ES unitsSiii4^. The difference 
in the intensities of these peaks indicates that the alloy 
is formed basically by CS tetrahedra and ES tetrahe- 
dra are found in a small quantity. The shoulder around 
237 cm -1 is associated with A\ and E modes of Se chains, 
and they are also seen in the c-Se (see Fig.QJ. The broad 
band at 255 cm -1 is related to A\ and E% modes of Se„ 
rings^&iiiiSi^. It is important to note that the band at 
around 165 cm" 1 , which is associated with Ge-Ge pairs 
vibrations in ethanlike units, is not seen in the spectrum 
of the alloy (as well as the other bands of c-Ge), indicat- 
ing that the quantity of Ge-Ge pairs is very low in the 
alloy. In addition, the bands related to Se-Se pairs are 
not seen in the spectra of the alloys produced by MQ 41 
or VE& techniques at this composition, but they exist in 
the alloy produced by MA, and their intensities suggest 
that their quantities may be relevant. These results in- 
dicate that the tetrahedral units are formed during the 
MA process and there is a preference for CS units. 
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FIG. 1: RS spectra of (a) MA-a-Ge 3 oSe 7 o, (b) c-Se and (c) 
c-Ge. 
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B. DSC measurement 

Figure shows the DSC measurement of the alloy, and 
it displays only two exothermic peaks at about 200°C and 
360° C, which are associated with crystallization temper- 
atures of MA-a-Ge 30 Se7o. There are no peaks associated 
with melting of c-Se, which should appear around 217°C, 
nor with glass transition or crystallization of o-Se, which 
should occur at about 45°C and 90°C, respectivery 5 * 4 ^. 
These data indicate that the Raman lines associated with 
Se chains and rings are not due to the presence of some 
quantity of unreacted a-Se or c-Se and reinforce the ex- 
istence of a high number of Se-Se pairs in the alloy, indi- 
cating that the local structure of MA-a-Ge 3 oSe 7 o proba- 
bly is different from that found in MQ, VE and even in 
MG samples, since now there are Se-Se pairs in a higher 
quantity. 
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FIG. 2: DSC measurement performed on MA-a-Ge3oSe7o. 



C. X-ray Diffraction and RMC Simulations 

Figure shows the experimental XRD S(K) (full line) 
for our alloy and the experimental ND S (K) (dashed 
line) given in Refs. Q or [3 for a-GeSe 2 . The FSDP is 
seen at around 1.1 A -1 . It is lower than those shown 
m Refs. and H indicating that the IRO in the alloy 
produced by MA is less pronounced than in the MQ- 
GeSe2 sample a-^- 4 . The FSDP is known to be much 
dependent on Ge-Ge and, to a lesser extent, on Ge-Se 
correlations 4 - 11 - 45 46 . Therefore these correlations have a 
different behavior in MA-a-Ge3oSe7o- 

S(K) was simulated using the RMC program available 
on the Interne1i2i. To perform the simulations we have 
considered a cubic cell with 16000 atoms (4800 Ge and 
11200 Se), 5 = 0.002, and a mean atomic number density 
po = 0.03868 ± 0.0005 atoms/A 3 . This value was found 
from the slope of the straight line (— Anp^r) fitting the 
initial part (until the first minimum) of the total G(r) 



functional. The minimum distances between atoms were 
fixed at the beginning of the simulations at tq c _q c — 
1.8 A, r Gc _ Sc = 1.8 A and r Sc . Sc = 1.75 A. The S RMC (K) 
obtained from the simulations (squares) is also shown 
in Fig. |3| and there is a very good agreement with the 
experimental S (K) . 
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FIG. 3: Experimental (full line) and simulated (squares) to- 
tal structure factor for MA-a-Ge3oSe7o together with the ND 
total structure factor given in Ref. (dashed line) (for a 
better comparison it was cut at K max = 8.5 A" 1 ). 

First hard sphere simulations without experimental 
data were carried out to avoid possible memory effects 
of the initial configurations in the results. Then un- 
constrained runs (i.e. when only hard sphere diameters 
and experimental data were used during the simulation) 
were carried out. These runs led to essentially identi- 
cal G RMC (r) functions, because of the proximity of the 
atomic numbers and scattering factors of Ge and Se. In 
the next series of simulations we tried to fit the exper- 
imental S(K) using the ND coordination numbers re- 
ported in Refs. Q or as coordination constraints. In 
this case we reached a poor agreement between the ex- 
perimental and the simulated structure factor, indicating 
that the local structure of MA-a-Ge3oSe7o is really dif- 
ferent from that found in the MQ-GeSe2 alloy, in agree- 
ment with Raman results. Thus, based on Raman results 
(small number of Ge-Ge pairs, presence of CS and ES-like 
units and a high number of Se-Se pairs) we used as start- 
ing coordination constraints the values TVgb-Gc = 0.2, 
A^Ge-Se = 3.7 and Ag c _s G = 1-0, which were then allowed 
to vary freely. The results obtained from the best simula- 
tion achieved are shown in Fig. [3] and they are discussed 
below. As a final test, we also tried to make simula- 
tions forcing the Se-Se coordination to be two, as it is 
in a-Se, but again the simulations did not show a good 
convergence. When this constraint is released, the Se-Se 
coordination number decreases until it reaches the num- 
ber we found for the previous case (As e -Sc = 1-25, see 
text below), in agreement with Raman and DSC results. 
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1. Pair Distribution Functions 

Figure 2] shows the G^ M (r) functions obtained from 
the RMC simulations of MA-a-Ge 30 Se7 compared to 
those found for MQ-GeSe^- There are important differ- 
ences among these functions and those obtained by ND 
measurements^ and also by MD simulationsi2*ii*i^*i&. 
First of all, we should note that the densities of the al- 
loys are very different (pmq-GcSc 2 = 0.0334 atoms/A -3 
— ), so the heights of the peaks in Gy (r) cannot be easily 
compared. The intensity of the first peak in G^g^ is 
higher, confirming their existence in a larger quantity in 
the first coordination shell, as indicated by the RS mea- 
surement. The two first peaks of the G™g c (r) function, 
which correspond to Ge-Ge first and second neighbors, 
show up around 2.33 and 3.83 A. The first peak occurs 
at a distance a little shorter than that found in the MQ 
sample 1 - 4 , but the second is displaced towards higher- 
r values by 0.26 A. Besides that, a minimum occurs at 
around 3.0 A. Remembering that the distance between 
two Ge atoms in adjacent ES and CS units are found 
at 3.02 A and 3.60 A, respectiveljiii*i^, it can be seen 
that the fraction of ES units in MA-a-Ge3oSe7o is low, in 
agreement with the results obtained by RS. If we consider 
that in c-GeSe2 each ES or CS Ge atom has, respectively, 
three and four nearest-neighbor Ge atoms and there are 
equal numbers of ES and CS units giving Nq c _q c = 3.5 
we find that, for a value of Ag c -Gc = 3.85 (see Table QJ, 
we should have about 85% of the tetrahedral units as CS 
units and 15% as ES units. Since the intensity of the 




FSDP in S(K) seems to be related to the quantity of ES 
tetrahedraiSiiii 4 * 1 ^, the low quantity of ES units in our 
alloy could explain the low intensity of the FSDP. 

The first peak of GQ^g c (r) function is located at 
2.35 A, as it happens to the MQ-GeSe2 samplesi^ii. 
This shell in our alloy is lower and broader than in MQ- 
GeSe2 but, due to the difference in densities, Ge-Se co- 
ordination numbers are almost the same in both alloys. 
The next peak appears at 3.84 A, and it is higher than 
that at 4.96 A. In the MQ-GeSe2 samples, there is a 
peak around 3.02 A which is not seen, or at least not re- 
solved, in MA-a-Ge3oSe7o, and there are peaks at 3.78 A 
(smaller) and 4.66 A (higher). In c-GeSe2 an ES Ge atom 
has Se neighbors in the range 4-6 < r < 5.3 A, and a CS 
Ge atom has Se neighbors in the range 4.0 < r < 4.8 AiS* 3 . 
Then, we have associated the peak at 3.84 A with CS 
units and that at 4.96 A with CS and ES units. 

The first peak of G^ s °(r) function is located at 2.33 A, 
and it corresponds to a coordination number ATgc-Sc = 
1.25, which is much higher than that obtained for the 
MQ-GeSe 2 samples (N Se -Se = 0.20)ii This suggests 
that some of the tetrahedral units are connected by Se 
"bridges" , forming sequences such as Ge-Se-Se-Se-Ge. As 
a consequence, Ge-Ge pairs should be found at higher 
distance values. This agrees with the results obtained 
from the G™Qe( r ) function. The next peak appears at 
3.75 A, which gives a ratio of Ge-Se:Se-Se distances of 
0.627. The valu e expected for ideal tetrahedral coordi- 
nation is y/3/8 = 0.612, indicating that the tetrahedral 
(CS or ES) units are distorted in our alloy. It is interest- 
ing to note that in MQ-GeSe2 samples no peaks are found 
from «4.5Ato«5.9A either considering ND results^ 
or MD simulations^*^. On the other hand, in MA-a- 
Ge 30 Se7o there are peaks at 4.95 and 5.75 A, and we 
believe these peaks are related to the distances between 
Se atoms in the "bridges" and Se atoms in tetrahedral 
units. The peaks at around 6.6 and 7.4 A, which are also 
seen in MQ-GeSe2 samples, can be associated, following 
Ref. ITTl to distances between Se atoms belonging to two 
adjacent tetrahedral units. 

The coordination numbers show in Table [I] were cal- 
culated considering the RDF?- MC (r) functions shown in 
Fig. The integrations were made using the following 
ranges: from 1.7 A to 2.95 A to the first peak and from 
2.95 A to 4.5 A to the second peak, for all RDF^ MC (r) 



functions. 
Tabled 



The interatomic distances are also shown in 



2. Partial Structure Factors 



FIG. 4: Ggo-GcM, Ggc-ScM and Gse-Se(r) functions ob- 
tained from RMC simulations (solid lines) and by ND (dashed 
lines, Ref. i). 



The partial Sjj MC (K) are shown in Fig. |H| together 
with the S^ D (K) found in Ref. H S^g e (K) has its 
first three peaks at about 1.1, 2.0 and 3.6 A -1 and two 
minima at 1.3 and 2.8 A -1 . Their positions agree rea- 
sonably well with those found for MQ-GeSe 2 studied by 
NDiii, but intensities are very different. The first peak is 
much lower in MA-a-Ge3oSe7o than it is in the MQ-GeSe2 
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MD simulationsiSiii was found only at 1.3 A 1 , where 
our Sq^g c (K) has a minimum. 



FIG. 5: RDFi^L(r-), RDFg&(r) and RDF^r) ob- 
tained from the RMC simulations. 




TABLE I: 
Ge3oSe7Q. 



Structural Parameters obtained for MA-a- 



RMC 



First Shell 


Second Shell 


Bond Type Gc-Gc Ge-Sc Sc-Go Se-Se 

N 0.26 3.50 1.75 1.25 
r (A) 2.33 2.35 2.35 2.33 


Gc-Gc Gc-Sc Sc-Gc Sc-Se 

3.85 7.4 3.7 9.7 
3.83 3.84 3.84 3.75 


MQ-GeSe 2 studied by NBii 


Bond Type Gc-Gc Ge-Sc Sc-Gc Sc-Se 

N 0.25 3.7 1.8 0.20 
r (A) 2.42 2.36 2.36 2.32 


Gc-Gc Gc-Se a Sc-Gc a Sc-Se 

3.2 - - 9.3 
3.57 - - 3.89 


2-GeSe 2 studied by NDi 


Bond Type Gc-Gc Ge-Sc Sc-Go Se-Se 

N 0.25 3.5 1.7 0.23 
r (A) 2.33 2.42 2.42 2.30 


Gc-Gc Gc-Se Sc-Go Sc-Se 

2.9 4.0 2.0 9.6 
3.59 4.15 4.15 3.75 



"These numbers are not given in Refs. 1 or 4. 



sample. This peak is associated with the FSDP in the 
S{K) shown in Fig. [3] and, since the FSDP is known to 
be strongly dependent on the Ge-Ge correlations and, to 
a lesser extent, on the Ge-Se correlationsAi^^, the low 
intensity FSDP experimentally observed in Fig. could 
be caused by weak Ge-Ge and Ge-Se correlations at its 
position. In addition, the heights of the second and third 
peaks are almost the same, and in MQ samples the height 
of the second peak is twice of that of the third peak. It is 
interesting to note that the first peak obtained from the 



FIG. 6: Sgc-Gc(K) , iSGe-Se(^0 and Sse-Sc(K) factors obtained 
from the RMC simulations (thick lines) and by ND (thin lines, 
Ref.i). 

S g£sL( k ) nas two maxima at 1.02 and 3.58 A" 1 and 
a minimum at 1.7 A -1 . Again, the first peak of MA- 
a-Ge3oSe7o is lower than that found in MQ-GeSe2 sam- 
ples, but its position is the same, and it is also associ- 
ated with the FSDP in S(K). The second peak and the 
first minimum in MQ samples are found at 3.5 A -1 and 
2.1 A -1 , respectively, indicating that in our alloy they 
are dislocated to higher (the maxima) and lower- if (the 
minimum) values. These facts can be explained by the 
important differences between the Gg^c?, (r) function and 
that of the MQ samples for r > 4 A. 

S^iP(K) has two peaks at around 2.0 and 3.6 A -1 , 
and there is a minimum at 2.8 A" 1 . At about 1.2 A -1 
there is a small peak associated with the FSDP in S(K). 
In MQ-GeSe2 samples the peaks are seen at 0.95, 2.05 
and 3.55 A -1 , and there is a minimum at 2.75 A -1 . 
These data indicate that Ss^i(K) is similar to that of 
MQ-GeSe2 samples, at least in the \ow-K region, con- 
cerning peak positions. However, we should note that 
their heights are different, in particular the intensity of 
the peak at 2.0 A -1 . 

It is important to compare our results considering the 
BT formalism. Figure □ shows the S^ C (K), S^ C (K) 
and S^ C (K) factors obtained using Eqs. □ and 
together with those factors found by ND and shown in 
Ref. H As expected, S^ C (K) resembles the XRD 
S(K) because the scattering lengths of Ge and Se are 
almost the same, and it is very similar to the S^(K) 
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found for MQ-GeSe2, except for the FSDP intensity at 
1.0 A -1 . Although there are differences in the peak inten- 
sities, S^ C (K) resembles that found for l-GeSe^i 12 ' 13 i 49 
and for MQ-GeSe2 — , including the sharp minimum at 

2.0 A -1 . Scc(K) also behaves like that obtained for l- 
GeSe 2 ^ and MQ-GeSe 2 -. In Z-GeSe 2 a sharp FSDP is 
clearly seen in Scc(K) at around 1.0 A" 1 , and this fact 
also occurs in MQ-GeSe2 samples. MD simulations of 
l-GeSe 2 ^£ and MQ-GeSe 2 — could not reproduce well 
this peak. In our case, a very weak FSDP can be seen at 

1.1 A- 1 MS M shoulder of the high peak at 2.0 A 1 . Re- 
membering that the FSDP in S(K) (see Fig. EI) has a low 
intensity and considering all the results discussed above, 
we should not expect that the FSDP in Scc(K) was as 
high and well defined as it is in I- or MQ-GeSe 2 samples. 
The IRO in MA-a-Ge3oSe7o is different mainly because of 
the introduction of Se-Se first neighbor pairs as "bridges" 
between the tetrahedral units, thus decreasing the num- 
ber of ES units and increasing the number of CS units. 
This affects the short and medium range order, which, 
in its turn, changes the concentration-concentration BT 
factor. 




K (A ) 

FIG. 7: Bathia- Thornton Skn(K), Shc(K) and Scc{K) fac- 
tors obtained from RMC simulations (solid lines) and by ND 
(dashed lines, Ref. |^). The arrow indicates the FSDP in 



3. Bond-Angle Distribution Functions 

By defining the partial bond-angle distribution func- 
tions @iji(cos9) where j is the atom in the corner we 
calculated the angular distribution of the bonds between 



first neighbor atoms. The six ©^/(cosfl) functions are 
shown in Fig. |SJ All these functions were calculated con- 
sidering as r max the position of the first minimum after 
the peak of the first shell (r max « 3.0 A). 
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FIG. 8: Qiji(cos8) functions obtained from RMC simula- 
tions. 



The ©Ge-Ge-Ge (cos 9) function is very noisy because of 
the very small number of Ge-Ge pairs in the first shell, 
but it shows a tendency for angles around 100°. The 
©Se-Se-Se (cos #) function has peaks at 55-61° and a broad 
distribution from 99 to 118°, with a maximum at 103°. 
The internal Se-Se-Se angles in perfect tetrahedra are 
found at 60°, the Se-Se-Se angle in trigonal Se is seen at 
103° and the angles in small Se chains and rings 4 ^ can 
be found at 90-116°. Thus, the @s e -Se-Se(cos#) func- 
tion indicates that the tetrahedra in MA-a-Ge3oSe7o are 
slightly distorted and Se chains and rings are formed, in 
agreement with RS results and with the previous analy- 
ses of the Gij(r) functions. The ©s e -Ge-Se(cos 9) function 
is very similar to the ©se-Se-Se(cos 9) function, showing 
peaks at 58° and 105°, which is close to the ideal tetra- 
hedral angle of 109°. 

The ©Ge-Se-Ge(cos 9) function peaks at about 58° and 
106°. The Ge-Se-Ge sequence in ES units has angles 
around 80° that are not seen in this function, reinforcing 
the small quantity of these units. On the other hand, 
this sequence in CS units has angles around 100°, and 
n-fold rings have angles ranging from 92° to 125° li 4 * 4 ^ 
The ©se-Gc-Gc(cos 9) and ©G e -Se-Se(cos 9) functions are 
similar to the others, having peaks at 58°, associated 
with threefold rings, and at 116° (Se-Ge-Ge) and 109° 
(Ge-Se-Se), which are related to tetrahedral angles and 
?i-fold rings. The ©^ ; (cos 9) functions above confirm that 
distorted tetrahedral units are formed in MA-a-Ge3oSe7o, 
with a clear preference for CS units. These units seem to 
be connected by Se-Se bridges, forming small chains and 
rings, as pointed out by the RS data. 
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V. CONCLUSION 

To summarize, we can conclude that the amorphous 
Ge 3 oSe 70 alloy can be produced by MA starting from 
the elemental powders of Ge and Se, but the structure 
of the alloy is different from that found in MQ-, VE- 
or MG-GeSe2 samples, making clear the importance of 
the preparation technique. Structural units similar to 
distorted CS and ES tetrahedra are formed, with a pref- 
erence for CS tetrahedra, as indicated by Raman spec- 
troscopy. These units seem to be connected by Se-Se 
bridges, as suggested by the high number of these pairs 
in the first shell, by the increase in the distance between 
Ge-Ge second neighbors, by the 0y((cos#) functions and 
also by RS vibrational data. These differences in the SRO 
affect the IRO of the alloy, and this causes the low in- 
tensity of the FSDP in S(K) of MA-a-Ge 3 oSe 7 o when 
compared to the MQ-GeSe2 allojii The low intensity 
of the FSDP in S (K) can be traced back to the partial 
Sg c -Gc(K) factor, whose FSDP is related to the ES units, 



which are found in a small quantity in the alloy. Sec (K) 
reflects these features and shows a very weak FSDP when 
compared to the factor found for l-GeSe-^^^ or MQ- 
GeSeA 

As a second remark, this study reinforces the rele- 
vance of using and combining RMC simulations with 
other techniques to model amorphous structures, since 
all features described above were obtained considering 
directly the experimental S(K) in the simulations. 
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